Mixed flux-equipartition solutions of a diffusion model of nonlinear cascades 
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We present a parametric study of a nonlinear diffusion equation which generalises Leith's model 
of a turbulent cascade to an arbitrary cascade having a single conserved quantity. There are three 
stationary regimes depending on whether the Kolmogorov exponent is greater than, less than or 
equal to the equilibrium exponent. In the first regime, the large scale spectrum scales with the 
Kolmogorov exponent. In the second regime, the large scale spectrum scales with the equilibrium 
exponent so the system appears to be at equilibrium at large scales. Furthermore, in this equilibrium- 
like regime, the amplitude of the large-scale spectrum depends on the small scale cut-off. This is 
interpreted as an analogue of cascade nonlocality. In the third regime, the equilibrium spectrum 
acquires a logarithmic correction. An exact analysis of the self-similar, non-stationary problem 
shows that time-evolving cascades have direct analogues of these three regimes. 

PACS numbers: 47.35.-i, 82.20.-w, 94.05. Lk 



Cascades are often observed in the non-equilibrium sta- 
tistical dynamics of interacting many-body systems in 
which microscopic interactions between degrees of free- 
dom are conservative and sources and sinks of the con- 
served quantity are widely separated. A famous example 
is the Richardson cascade in high Reynolds number hy- 
drodynamic turbulence. There the nonlinear terms in the 
Navier-Stokes equation conserve the fluid kinetic energy 
and energy injection (by stirring for example) and energy 
dissipation (by viscosity) are widely separated in scale 
or wavenumber. See [l[ and the references therein. In 
recognition of the historical importance of this example, 
and for brevity, we shall always speak in this article of 
a cascade as the process whereby nonlinear interactions 
conservatively transport "energy" , E, in "wavenumber" 
space, k. We acknowledge that cascades occur in many 
other contexts in which the conserved quantity is not nec- 
essarily the energy and transport may occur in a space 
other than the space of wave-numbers. Some examples 
include wave turbulence 2] , cluster-aggregation Q , non- 
linear diffusion [4] and Bose-Einstein condensation 0, Q . 

The description of cascades based on the underlying 
dynamical equations typically leads to Boltzmann-like 
kinetic equations obtained by moment closures which 
may be phenomenological (as is common with hydrody- 
namic examples) or which may be asymptotically exact 
(as with weak wave turbulence Transport in such ki- 
netic equations usually involves a nonlinear integral col- 
lision operator. This makes their detailed analysis diffi- 
cult. One way around this difficulty, originally proposed 
by Leith Q, is to phenomenologically replace the inte- 
gral collision operator with a more analytically tractable 
nonlinear differential operator in such a way as to pre- 
serve the scaling properties of the original problem and, 
in particular, the scalings of the stationary Kolmogorov 
and equilibrium solutions. Such models are referred to 



as differential approximation models. These phenomeno- 
logical models allow the cascade dynamics to be qual- 
itatively explored with relative ease. For that reason, 
they have become an active field of research in their own 
right. The Leith model continues to be of interest in 
turbulence [§], while similar models have been used to 
study two-dimensional turbulence l3, Uj, wave turbu- 




13|, kelvin waves on vortex lines in a super- 



|14| . the Boltzmann equation for a hard-sphere gas 
and optical turbulence [l^. They are often used 
as a heuristic way of establishing the direction of the 



Kolmogorov cascade [15|, [16[ , an issue which has caused 
controversy in some contexts. As we shall see below, 
care must be taken to correctly interpret the predictions 
about the cascade direction made using these models. 

A disadvantage of differential approximation models is 
that, by construction, they model the nonlinear trans- 
port as a process which is local in scale. It is well known, 
however, that the integral character of the original colli- 
sion integral may lead to cascade dynamics which are 
nonlocal in scale. That is to say the energy transfer 
through a given wavenumber is dominated by interac- 
tions with smallest or largest wavenumber in the system 
rather than with nearby wavenumbers. The phenomeno- 
logical nature of differential approximation models opens 
the door to the uncomfortable possibility of attempting 
to model a nonlocal cascade with a local operator. 

Consider the following generalisation of Leith's model: 



dE_ dJ_ 
'dt ^ ~'dk' 
where J is the energy flux which is modeled as 

a 



J{k) 



dk 



(r^ E). 



(1) 



(2) 



The sign of the flux is chosen such that the flux is posi- 
tive when the energy flows to the right in fc-space. This 
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is easily verified by integrating Eq. ([T]) from to K (as- 
suming tliat J(0) = 0) and asking whether energy is en- 
tering or leaving the interval [0,-ftr]. Eq. (HI) has three 
adjustable parameters, m, xk and xt- The parameter 
771, which we take to be greater than 1, is the order of 
the nonlinear interaction responsible for the transport of 
energy. The parameters xk and xt are, as we shall see 
below are the exponents of the stationary Kolmogorov 
and thermodynamic equilibrium states, which we take to 
be independent adjustable parameters in order to per- 
form a parameteric study of the properties of the gener- 
alised Leith model. Note that the original Leith model 
is recovered by setting m = 3/2, xk=5/3 and xt = —2. 
Eq. Il]) is appropriate for modeling a system with a single 
conserved quantity, and thus a single cascade. We should 
remark at this point that Eq. ^ is also of considerable 
independent interest outside of its utility as a heuristic 
model of turbulent cascades. For various values of the 
parameters, it appears as a model of flow in a porous 
medium viscous gravity currents 17 1, transport of 
density fluctuations in a magnetised plasma Il8| and the 
spreading of surfactant on a liquid interface 19*. 

The general stationary solution of Eq. ([T]) involves two 
constants which we call J and T: 



rj-rn 



J 



XK - XT 



(3) 



There arc two stationary solutions which are pure power 
laws. The first, having J = is 



E{k) ^Tk- 



(4) 



It corresponds to the thermodynamic equilibrium solu- 
tion since the flux, Eq. ([2]), vanishes on this solution. 
For this reason, we refer to T as the temperature even 
though, strictly, the thermodynamic temperature is only 
defined at equilibrium. The second, having T = 0, is 



J 



XK - XT 



k- 



(5) 



It corresponds to the Kolmogorov solution since the flux, 
Eq. ([2), is constant and equal to J on this solution. 

From Eq. ([S|) one can see that the flux, J, carried by the 
Kolmogorov spectrum must be positive when xk > xt 
corresponding to energy transfer to the right in fc-space 
and negative when xk < xt corresponding to energy 
transfer to the left in fc-space. Energy transfer to the 
left is inconsistent with the energy injection occuring at 
small fc and the energy dissipation occuring at small fc 
- that is the flux is in the "wrong direction" to connect 
the source and sink. The identification of situations in 
which this happens is one of the popular applications 
of differential approximation models. The issue is not 
entirely theoretical and occurs in reality for energy and 
particle cascades in the Boltzmann equation (20| and for 



the inverse cascade in two-dimensional optical turbulence 
(lif . It is generally agreed that this means that the Kol- 
mogorov spectrum is not physically realisable although 
there is less consensus about what takes its place. We 
address this issue clearly and unambiguously below, at 
least in the context of the generalised Leith model. 

Let us return to the general stationary state, Eq. ([31 
which has finite J and T. There are three cases: 
1 Kolmogorov-like regime, xk > xt: In this regime, 
we have a regular Kolmogorov cascade at large scales: 



J 



^xk — Xt , 
which is thermalised at small scales: 



as fc 0, (6) 



E{k) ^Tk- 



as fc — ^ oo. 



(7) 



Such states are sometimes called "warm" cascades [21[ 
since they have a nonzero temperature parameter. They 
are relevant for the description of the statistical dynamics 
of the truncated Euler equations for example [22l|. 
2 Equilibrium-like regime, xk < xt: In this regime, 
the cascade has a completely different character. It ap- 
pears to be at equilibrium at large scales, despite carrying 
a constant flux: 



E{k) ^Tk- 



as fc — 0. 



(8) 



Note that J can be positive in this regime. Eq. ([3]) can 
therefore describe a cascade with the "correct" direction 
provided we allow a flnite value of T . The reason is that 
the term describing the flux is subleading. In contrast to 
the case xk > xt, the cascade has an intrinsic cut-off at 
which the spectrum vanishes given by 



^ XT 



XK 



J 



T 



(9) 



It is, perhaps, more natural to consider the temperature, 
T, as a function of the cut-off, fc*, which may be imposed 
for example by the dissipation scale. In this case, the 
stationary state can be written as 



J 



E{k) = k- 

XK - XT 

(10) 

The amplitude of the spectrum at large scales depends on 
the small scale cut-off in this regime. This is the analogue 
of non-locality for the differential approximation model. 
Intriguingly, it was shown in [23[ that the isotropic 3- 
wave kinetic equation is always nonlocal when xk < xt- 
Thus not only does Eq. (fTO|) provide us with an analogue 
of non-locality for Eq. ([T]) but it occurs for the correct pa- 
rameter regime. A relationship between the temperature 
and the small-scale cut-off has recently been proposed 
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and partially observed numerically in the context of the 
classical Boltzmann equation fl5|. 

3 Degenerate regime, xk = xx: When the thermody- 
namic and Kolmogorov exponents coincide, the spectrum 
is: 



E{k) = (Jm)- k- 



(11) 



so that the system appears to be at equilibrium with a 
logarithmic correction. This can be seen by direct inte- 
gration of the stationary version of Eq. ((TJ. It is more 
informative to obtain this formula by Taylor expanding 
Eq. ()10p in small values of xk — xt and taking the limit 
Xk Xt- By doing this, one sees clearly that the log- 
arithm is the remnant of the cut-off dependence of the 
equilibrium-like regime as one enters the cut-off indepen- 
dent Kolmogorov-like regime. While this case is a special 
point in the parameter space, it does occur in practice as 
for example in the direct cascade in 3-D NLS turbulence 
16| and in elastic wave turbulence in a vibrating plate 

Let us now consider non-stationary cascades relevant 
to situations in which we do not inject energy at large 
scales but rather consider the evolution of a lump of en- 
ergy which is initially concentrated at large scales. In this 
case, there is no stationary spectrum and the evolution 
is described by a self-similar function of k and t: 



E{k,t) ^ s{tr F{0 



^here ^=^, 



(12) 



where s(t) is a typical wavenumber which grows in time 
as the spectrum spreads in fc-space. It is well-known (see 
for example [1^ and the references therein) that this self- 
similarity ansatz applied to nonlinear diffusion equations 
like Eq. ^ leads to weak solutions describing propagat- 
ing fronts which are positive on an expanding compact 
interval, [0,fc*(t)], and zero elsewhere. It is convenient 
therefore to take the characteristic scale, s(t), to be the 
right boundary of the support of the solution correspond- 
ing to the front tip. We show now that the same three 
regimes identified above for the stationary case have di- 
rect analogues in the non-stationary case. Substituting 
Eq. (fT^ into Eq. (H]) we obtain the scaling equations: 



aF 



ds 
di 

dF 



d 



-xt + 1 pr> 



d^ 



(13) 
.(14) 



of energy, J^^"'"*^ E{k,t) dk — 1 
^ F(^) = 1, from which we 



From conservation 

Eq. ([HI) leads to s 
conclude that a = — 1. To avoid the complications 26 1 
associated with self-similarity of the second kind, we shall 
assume from this point on that xk < 1- Thus the results 
cited below are applicable to infinite capacity cascades 
only. When a = — 1, the left hand side of Eq. ((Ti)) is 



an exact differential and can be integrated explicitly to 
obtain the scaling function in closed form (assuming that 
1 — mxK + {ni — ^)xt ^ 0): 



F{0 = 



I ^ — XT 



xt-xns 



if < e < 1 
otherwise 

(15) 

where we have introduced, for convenience, the nonsta- 
tionary exponent 



xns 



m Xk — 1 
TO — 1 



(16) 



This solution is the analogue for the generalised Leith 
model of the front solutions of the porous medium equa- 
tion originally obtained by Pattle 27]. As before, exam- 
ining Eq. ()15|) shows that there are 3 regimes: 
1 Nonstationary Kolmogorov-like regime, 
Xns > In this regime, we have a non-stationary 
cascade at large scales with the Kolmogorov exponent, 
Xk, replaced by xns'- 



1 



Xns - XT 



as ^ ^ 0. (17) 



2 Nonstationary equilibrium-like regime, 

Xns < Xx: In this regime the cascade appears to 
be at equilibrium at large scales, but with a temperature 
which decays in time due to Eq. (IT2t : 



1 



XT - XNS 



as C ^ 0. (18) 



3 Nonstationary degenerate regime, xk = xx: 

When the nonstationary and equilibrium exponents co- 
incide, we obtain the non-stationary analogue of the log- 
arithmic correction to the equilibrium scaling discussed 
above for the stationary case: 




lOf 



if < C < 1 
otherwise 

(19) 

Note that, as in the stationary case, the flux is posi- 
tive and to the right in all three regimes. We could also 
consider nonstationary cascades with a source of energy, 
in which case, the exponent a would no longer be equal 
to —1 and we would lose the exact differential on the 
left hand side of Eq. ([M)) which allowed us to solve the 
problem explicitly. In this case, however, the asymp- 
totics of the scaling function, F{£), can be obtained using 



the phase plane methods developed in [17|. This rather 
technical analysis will be presented elsewhere. In order 
to connect the results presented here on the generalised 
Leith model back to the integral collision operators which 
we purport to model, we remark that the nonstationary 
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regimes 1 and 3 have already been explored in consid- 
erable detail analytically and numerically in the context 
of the isotropic three- wave kinetic equation [2^ and con- 
form to the general behaviour outlined here. 

To conclude, we have performed a complete paramet- 
ric study of a nonlinear diffusion model of a turbulent 
cascade with a single conserved quantity which gener- 
alises Leith's original model of the energy cascade in 3 
dimensional hydrodynamic turbulence. Both stationary 
and non-stationary cascades can be described by simple 
analytic solutions of the model. We showed that there 
are three regimes depending on whether the Kolmogorov 
exponent is greater than, less than or equal to the equilib- 
rium exponent. In the Kolmogorov-like regime, the equi- 
Ibrium behaviour is a small correction to the finite fiux 
spectrum at large scales. Large scales are independent 
of the small scale cut-off. In the equilibrium-like regime, 
the finite flux behaviour is a small correction to equilib- 
rium spectrum at large scales. The amplitude of the large 
scale spectrum is a diverging function of the small scale 
cut-off. This latter fact means that even differential ap- 
proximation models can mimic some aspects of cascade 
nonlocality. In the degenerate regime, both finite fiux 
and equilibrium behaviours are equally important leading 
to a logarithmic correction to the equilibrium spectrum. 
The question of cascade direction is completely clear in 
this model. By allowing a finite T, the flux is always pos- 
itive and in the "correct" direction. The Kolmogorov-like 
and degenerate regimes are already well known but the 
equilibrium-like regime has not been appreciated previ- 
ously and should now be sought in kinetic equations using 
the full collision integral. Finally, we remark that many 
of the interesting physical applications of differential ap- 
proximation models in which issues of cascade direction 
and degeneracy of exponents arise have two conserved 
quantities. This considerably complicates things because 
the corresponding differential equation is fourth order. 
We hope that the comprehensive description of a single 
conservation law presented here will help to clarify the 
issues arising in more complicated examples. 



ACKNOWLEDGEMENTS 



We acknowledge useful discussions with D. Proment 
and S. Nazarenko. 



* lconnaughtonc@g mail.coml 
t rm837@york.ac.uk] 

[1] U. Frisch, Turbulence: the legacy of A. N. Kolmogorov 
(Cambridge University Press, Cambridge, 1995). 

[2] V. Zakharov, V. Lvov, and G. Falkovich, Kolmogorov 
Spectra of Turbulence (Springer- Verlag, Berlin, 1992). 

[3] F. Leyvraz, Phys. Reports 383, 95 (Aug. 2003). 

[4] J. L. Vazquez, The Porous Medium Equation - Mathe- 
matical Theory (Oxford University Press, 2006). 

[5] R. Lacaze, P. Lallemand, Y. Pomeau, and S. Rica, Phys- 
ica D 152-153, 779 (2001). 

[6] C. Connaughton and Y. Pomeau, C. R. Phys. 5, 91 
(2004). 

[7] A. C. Newell and B. Rumpf, Ann. Rev. Fluid Mech. 43, 
59 (2011). 

[8] C. E. Leith, Phys. Fluids 10, 1409 (1967). 
[9] T. Clark, R. Rubinstein, and J. Weinstock, J. Turbulence 
10, 1 (2010). 

[10] C. E. Leith, Phys. Fluids 11, 671 (1968). 

[11] V. Lvov and S. Nazarenko, JETP Lett. 83, 541 (2006). 

[12] S. Hasselmann, K. Hasselmann, J. H. AUender, and T. P. 
Barnett, J. Phys. Ocean. 15, 1378 (1985). 

[13] C. Connaughton, A. Newell, and Y. Pomeau, Physica D 
184, 64 (2003). 

[14] S. Nazarenko, JETP Lett. 83, 198 (2006). 

[15] D. Proment, M. Onorato, P. Asinari, and N. Nazarenko, 
"Nonequilibrium steady solutions of the Boltzmann equa- 
tion," (2010), arXiv:1101.4137vl [nlin.CD]. 

[16] S. Dyachenko, A. Newell, A. Pushkarev, and V. Za- 
kharov, Physica D 57, 96 (1992). 

[17] J. Gratton and F. Minotti, J. Fluid Mech. 210, 155 
(1990). 

[18] V. Naulin, A. H. Nielsen, and J. J. Rasmussen, Phys. 

Plasmas 12, 122306 (2005). 
[19] O. E. Jensen, J. Fluid Mech. 293, 349 (1995). 
[20] A. V. Kats, V. M. Kontorovich, S. S. Moiseev, and V. E. 

Novikov, ZhETF Pis ma Redaktsiiu 21, 13 (1975). 
[21] C. Connaughton and S. Nazarenko, Phys. Rev. Lett. 92, 

044501 (2004). 

[22] C. Cichowlas, P. Bonati, F. Debbasch, and M. Brachet, 

Phys. Rev. Lett. 95, 264502 (2005). 
[23] C. Connaughton, Physica D 238, 2282 (2009). 
[24] G. Diiring, C. Josserand, and S. Rica, Phys. Rev. Lett. 

97, 025503 (2006). 
[25] A. A. Lacey, J. R. Ockendon, and A. B. Tayler, SIAM J. 

Appl. Math. 42, 1252 (1982). 
[26] C. Connaughton and A. C. Newell, Phys. Rev. E 81, 

036303 (2010). 

[27] R. E. Pattle, Quart. J. Mech. Appl. Math. 12, 407 (1959). 
[28] C. Connaughton and P. Krapivsky, Phys. Rev. E 81, 
035303(R) (2010). 



